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Abstract 

We make a careful study about the nonrelativistic reduction of one-meson-exchange models 
for the nonmesonic weak hypernuclear decay. Starting from a widely accepted effective 
coupling Hamiltonian involving the exchange of the complete pseudoscalar and vector meson 
octets (it, n, K, p, ui, K*), the strangeness-changing weak AiV — > NN transition potential 
is derived, including two effects that have been systematically omitted in the literature, 
or, at best, only partly considered. These are the kinematical effects due to the difference 
between the lambda and nucleon masses, and the first-order nonlocality corrections, i.e., 
those involving up to first-order differential operators. Our analysis clearly shows that the 
main kinematical effect on the local contributions is the reduction of the effective pion mass. 
The kinematical effect on the nonlocal contributions is more complicated, since it activates 
several new terms that would otherwise remain dormant. Numerical results for and ^He 
are presented and they show that the combined kinematical plus nonlocal corrections have an 
appreciable influence on the partial decay rates. However, this is somewhat diminished in the 
main decay observables: the total nonmesonic rate, T nm , the neutron-to-proton branching 
ratio, T n /T p , and the asymmetry parameter, q,a- The latter two still cannot be reconciled 
with the available experimental data. The existing theoretical predictions for the sign of a\ 
in ^He are confirmed. 
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1 Introduction 



The free decay of a A hyperon occurs almost exclusively through the mesonic mode, A — ► irN, 
with the nucleon emerging with a momentum of about 100 MeV/c. Inside nuclear matter (pp ~ 
270 MeV/c) this mode is Pauli blocked, and, for all but the lightest A hypernuclei (A > 5), 
the weak decay is dominated by the nonmesonic channel, AN — > NN, which liberates enough 
kinetic energy to put the two emitted nucleons above the Fermi surface. In the absence of 
stable hyperon beams, these nonmesonic decays offer the only way available to investigate the 
strangeness-changing weak interaction between hadrons. (For reviews on hypernuclear decay, 
see Refs. [l]— |SJ.) 

The simplest model for this process is the exchange of a virtual pion and in fact this can 
reproduce reasonably well the total (nonmesonic) decay rate, T nm = T n + T p , but fails badly 
for other observables like the ratio of neutron-induced (An — * nn) to proton-induced (Ap — > np) 
transitions, T n /T p , and the asymmetry parameter a\. The deficiency of this model is attributed 
to effects of short range physics, which should be quite important in view of the large momentum 
transfers involved (~ 400 MeV/c). Although there have been some attempts to account for 
this fact by making use of quark models to compute the shortest range part of the transition 
potential most of the theoretical work opted for the addition of other, heavier mesons in 

the exchange process ^01 _ |2S]- None of these models gives fully satisfactory results. Inclusion 
of correlated two-pion exchange has not been completely successful either [241 125| . Nor have the 
addition of uncorrelated two-pion exchange, two-nucleon induced transitions or medium effects, 
treated within the nonrelativistic IMPISH or relativistic (32] propagator approaches, been of 
much help. 

Here, we concentrate on the line of one-meson-exchange (OME) models (l]-j25j. We do not 
explicitly discuss hybrid models, i.e., those involving quark degrees of freedom [3] - !II]- However, 
much of the theoretical developments we present could be generalized to include them. Also 
two-pion exchange |241 125j could be brought into our general framework. The main ingredients 
of OME models are the effective baryon-baryon-meson weak and strong Hamiltonians. These are 
constructed in the language of relativistic field theory, but in almost all calculations (exceptions 
are Refs. |14|-|16j) one has proceeded to make a nonrelativistic reduction for the extraction 
of the transition potential. This often involves some further approximations, like neglecting 
the nonlocality in this potential, and balancing by hand the distorted kinematics in the OME 
Feynman amplitudes, resulting from the difference between initial and final baryon masses. This 
not only alters the different terms in the transition potential, but also eliminates several of them. 
Our main purpose here is to assess the relative importance of these effects. 

The paper is organized as follows. Most of the formalism is developed in Section [21 In 
Subsection l2.il we explain the construction of the nonrelativistic transition potential, taking due 
care of the kinematics and including nonlocal terms, and, in Subsection l2.21 we give a motivation 
for not neglecting a priori the lambda-nucleon mass difference. In Subsections 12.31 and 12.41 
the explicit expressions for the local and first-order nonlocal contributions to the transition 
potential due to the exchange of a it, p, K or K* meson are derived and commented. The 
ones corresponding to the r/ and u exchanges can be easily obtained by analogy with those of 
the 7r and p, respectively, thus allowing the inclusion of the full pseudoscalar and vector meson 
octets, as deemed necessary by the present day consensus. In Subsection 12.51 we describe how 
to take the finite size effects into account by means of form factors. Our numerical results 
are reported in detail and discussed in Section [31 The phenomenological way to include short 
range correlations is presented in Subsection 13.11 together with the main expressions for the 
calculation of the transition rates in the extreme particle-hole model of Ref. |21j . and all this 
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is applied to the decay of ^C. In Subsection E2J we do the same for |He and compute also 
the asymmetry parameter. Finally, Section ^ summarizes our main conclusions. Some useful 
formulas are collected in Appendix EI and a specific point relating to our phase conventions is 
discussed in Appendix iBl 



2 OME transition potential 

2.1 General discussion 







/N 




Pi, s 'i\ 






x i (a /WW 


VW\(T) x 2 




Pi, Si/ 


\P2,S 2 




tJ 


w 



Figure 1: OME Feynman amplitude in coordinate space. 



The transition rate for the nonmesonic weak decay of a hypernucleus in its ground state |J), 
having energy Ej, to a residual nucleus in any of the allowed final states (F\, having energies 
Ep, and two outgoing nucleons is given by Fermi's golden rule, 



T = 2, £ 



s' s'„ F 



d 3 p[ d 3 p' 2 



(2tt) 3 (2tt) j 



d(E[ + E' 2 + E F - Ei) (M p' 2 s' 2 , F\V\I) 



(1) 



To construct the transition potential V in one-meson-exchange models, one starts from the free 
space Feynman amplitude depicted in Fig. ^ where x = (t, x) denotes space-time coordinates, 
p, momentum, and s, spin and, eventually, other intrinsic quantum numbers (such as isospin). 
In the remainder of this subsection we will consider a general situation, i.e., without specifying 
which baryons are propagating in each of the four legs, or which meson is being exchanged, 
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or yet the exact nature of the couplings at the two vertices. This will be particularized to 
A-hypernuclear decay in the subsections that follow. 

Vertices a and b correspond to coupling Hamiltonians of the general form (c = a or b) 

n c (x)=gMx)[T c (d)^(x)}i;(x), (2) 

where tp and <f> stand for the baryon and meson fields, respectively, g c is a coupling constant 
and r c may contain differential operators, in which case they are understood to be acting on 
the boson field only. The Feynman rules give 

T = (2tt) 4 5{E[ + E' 2 -E l - E 2 ) 5(p[ + p' 2 - Pl -p 2 )F, (3) 



where E t = ^Mf+pf (i = 1,2) for the incoming baryons, having masses Mj, and similarly 
(primed quantities) for the outgoing ones, and F is the Feynman amplitude in momentum space. 
Choosing the CM frame, 

- Pi = P2 = P> -Pi = P2 = P', ( 4 ) 

this can be put in the form 

iF(p', p; si4siS2) = x\ x\ V(p',p) Xai Xs 2 (5) 

with 

V(p',p) = u[(-p')u 2 ( P ') V(q) Ul (-p)u 2 (p), (6) 
where an d X ' ^i(Pi) (* = 1)2) are the momentum eigenspinors and their conjugates 

i 

for the incoming and outgoing baryons, respectively, and, denoting the meson propagator by T>, 

- iV(q) = [-ig a T a (iq)} iV{q) [-ig b T b (-iq)}. (7) 
We have introduced the 4- momentum transfer in the CM frame, q = (uj,q), with 

u = -{E 1 -E[ + E' 2 -E 2 ), q = p'-p. (8) 

Notice that we have directed q from vertex a to vertex b. 

The nonrelativistic transition potential V is given by the identification 

(p'\V\p) = V(p',p), (9) 

where an expansion up to quadratic terms in momentum/mass is implied. To this end it is 
convenient to change the momentum variables to q, defined in Eq. (jH]), and 

„ m p' + m'p , . 

Q = r- , , io 

m + m' 

where m and m' are the initial and final reduced masses, 

111 111 

m~W l + W 2 7ri~M[ + M 2 ~' ^ ' 

In this transformation, the following relations hold: 

P' 2 P 2 1 / 1 1 \ ^ Q 2 



2m' 2m 2 \m' mj 2(m + m') 
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and 



where 



p' ■ r' — p ■ r = Q ■ (r' — r) + q 



r 2 — T\ and 



mV + mr 
m + m' 



(13) 
(14) 



are, respectively, the initial and final relative coordinates. 

The contribution of any given meson, i, to Eq. (JJ3) has the general form 



Vi(p',p) 



Vi(p',p) 



q 2 + ri 



(15) 



where in is the meson mass. The nonrelativistic expansion of the numerator in Eq. f)15j) poses 
no problem, but the denominator does not truly have such an expansion. 1 Therefore, it needs a 
special treatment. Recall that, strictly speaking, the Feynman amplitude F in Eq. (J5J) is defined 
only for energy-conserving transitions, i.e., for E\ + Ei = E[ + E' 2 , and in this case one has, in 
the nonrelativistic approximation, 



PL 

2m 



P 



2m' 



M[ + M' 2 - Mi- M 2 . 



(16) 



This relation together with Eq. (J12I) allow us to write, again in the nonrelativistic approximation, 

g 2 Q 2 



U) 







2M„ 2Mr 



where we have introduced the kinematical masses Mq, M q and Mq, given by 

'Mi 



(17) 



Mi 







1 

1 + ^ 



M 2 Mi 
+ 



M 2 



2 VMi +M 2 M{ + M 2 



1 

+ 2 



1 _ 1 {Mt 



2 VMi + M 2 
1 (M x - M 2 M{ 



M 2 M\ 
+ 



M 2 



M{ 



M' 2 



(Mi - M() 
(M 2 -M 2 ), 



M„ 



M 2 



4 \Mi + M 2 M{ + My m + m' ' 



A/, 



M[ 



M< 2 



1 /Mi - M 2 
4 VMi + M 2 M{ + M> 2 



1 1 



/7? 



(18) 



It is clear that, whenever the absolute values of the differences of baryon masses are much smaller 
than the corresponding sums, as is the case, e.g., for hypernuclear decay, one can take advantage 
of the inequalities \Mq/M q \ <^ 1 and \Mq/Mq\ <C 1 to write the following approximation 



1 



1 



Mq 
Mq 



1 + Mi q' + m 



(1 + 



(q 2 + V 



2\2 



(19) 



where we have introduced the effective meson mass 



\ 1 1 Mo 

\ Mq 

x The nonrelativistic approximation is for the baryon dynamics. 



Ml 



(20) 



5 



The net result is that the nonrelativistic approximation, as defined here, will reduce V(p',p) 
to a quadratic polynomial in Q, 

V(p', p) - V® (q) + V W (g) • Q + Q • V( 2 ) (g) • Q, (21) 

whose coefficients , and V^ 2 ** are themselves, excluding the denominators that come 

from the meson propagators, polynomials of degree 2, 1 and 0, respectively, in q. This Q 
dependence translates, in the coordinate representation, into an expansion in the nonlocality of 
the transition potential. To see this, we make use of Eq. © to write 2 

d 3 p' f d 3 p 



j j^e*- r 'e-^ r V{p',p). (22) 



(2tt) 3 J (2tt) 

Changing the integration variables to (q,Q), making use of Eq. (|13|) and truncating, for sim- 
plicity, the quadratic polynomial in Eq. (|21ll at the linear term, this gives 

(r'\V\r) * [y(°)(r') - iV^ • V'] 5{r' - r), (23) 

where 

d*q 



and 



V {0) (r)= /^e^7(°)(g) (24) 

V«(r) = |-^ e *'V«(g). (25) 

This means that for any state of relative motion, VP, 

(r|y|^) = y(r)*(r), (26) 

with the transition potential in coordinate space, V(r), given, as an operator in wave- function 
space, by 

V(r) =v(°)(r) + ^)( r ) (27) 
where (r) is the local potential, and the differential operator 

y(i) (r) = _j m ( v . yW^)) _ ivWfr) • V, (28) 

its first-order nonlocality correction 3 . For our purposes here it will be sufficient to stop at this 
order, and we will not consider the second-order corrections, that would come from the last term 
in Eq. (ISTD . 

Before closing this subsection, let us make a brief comment on our treatment of the uj 2 term 
in the meson propagator in Eq. ()15j) . more specifically, on the expression we use for the time 
component, to, of the 4-momentum transfer in Eq. (JSJ). We are extracting the transition potential 
from the Feynman diagram in Fig. in which the baryons are on their mass shells. However, as 

2 To do this we have to assume that Eqs. ® and (12 1(1 hold for p and p' unrestricted, although V(p' ,p) was 
extracted from the Feynman amplitude F in Eqs. © and J^J, which, being related to a T-matrix, is unambiguous 
only for energy-conserving transitions. 

3 Notice that, despite its name, this has a local piece, namely, the first term in Eq. (128(1 . where (V • V^(V)) = 
divV (1) (r). 
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already alluded to in footnote El we need to extend these OME amplitudes to the off-shell region, 
for which there is no unique procedure. In the meson-exchange theory for the strong NN force, 
this is done by treating the two interacting particles through the Bethe-Salpeter equation 33 , 
and this ambiguity appears in the choice of which one of its various tridimensional reductions 
to use. This issue, which is related to meson retardation effects, has been much discussed in 
the past jSl] - !!^' but remains unsettled. We have followed the general philosophy of Machleidt 
and collaborators [HE] - !!!]]; to the effect that, for the case of similar masses, the best choice is to 
use tridimensional reductions that treat the two particles symmetrically, like the Blankenbecler- 
Sugar [32] or the Thompson equations. One, then, puts the two interacting particles equally 
off-shell in the CM frame, and fixes the time-components of the relative 4-momenta in the 
initial and final two-particle propagators at the values po = \{E 2 — E\) and p' = \{E' 2 — E[), 
respectively. (See, for instance, Eq. (2.28) in Ref. [HZ], remembering our convention in Eq. (J3J)-) 
This leads directly to our expression for uj in Eq. (J8J). Notice that, for strictly equal masses, this 
gives u = 0, being, therefore, equivalent, in this case, to the instantaneous approximation. Our 
choice for uj in Eq. @ leads, in the nonrelativistic approximation, to Eq. (|17|) and, consequently, 
to the expansion of the meson propag ators inEq. (JEJ). We are confident that this is appropriate 
for processes not too far off the energy-shell. In other situations that might occur, for instance, 
in a fully microscopic treatment of short-range correlations, this point should be reexamined. 



2.2 Kinematical effects 

In computing the OME Feynman amplitudes contributing to the strong NN force it is standard 
practice |39| to avoid the kinematical complications due to the difference between the neutron 
and proton masses by setting 

M n , M p -> M = (M n + M p )/2, (29) 

which can be justified by the small value of the ratio 

M n - M v 

" p = 0.0014. (30) 

M K ' 

The analogous practice is followed in the calculation of the transition potentials for the weak 
decay of A-hypernuclei ^H]. In this case, however, one equally sidesteps the lambda-nucleon 
mass difference, by setting, at the vetex where the A decays, 

M A ,M -> M = (M A + M)/2, (31) 

despite the fact that the corresponding ratio, 

M\ — M 



M 



0.17, (32) 



is nowhere as small. 

Undoubtedly, this approximation considerably simplifies the calculations. However, in view 
of the nonnegligible value of the ratio ([32 j). it seems appropriate to investigate the effects of 
the latter approximation. To this end, we examine below, for each meson in the pseudoscalar 
and vector octets, the expression for the nonrelativistic OME transition potential obtained by 
accepting the approximation in Ea. (|29|) . but not that in Ea. ([31|) . This gives for the kinematical 
masses (jlHj) 

M = \ (3M A + M) = 92.18 MeV , 



4 VM A + M 
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1 / M A 

2 \3M/ 

1 /M A - M 



M 



M q 2 V3M A + M 



4- = 2.196 x 10" 

M 



MeV~ 



8.800 x 10~ 5 MeV 



r-l 



(33) 



1 1 / Ma — M\ / 3Ma + M N 
Mq ~ 4 VMa + mJ V M A M 

where we have used @Ij M = 938.92 MeV and M A = 1115.68 MeV . The approximation (|5T|) 
would have set Mq, 1/M g and 1/Mq to zero. 



2.3 Contributions of nonstrange mesons 

For the nonstrange mesons, we have, acting respectively at the vertices a and b in Fig. ^ weak 
(W) and strong (S) coupling Hamiltonians that we take to be the same as those in Ref. |18j . 
For the pion they are 



Ti 



w 

ANn 



iG F nl ipN(A n + B n j 5 ) <f> 



V'A, 



(34) 



nf 



LNNn - WNN n^NlS W • T1p N , (35) 

where G F p% = 2.21 x 10 -7 is the Fermi weak coupling-constant, A n and B n are fitted to the 
free A decay, and gNNn is taken from OME models for the strong NN force. The isospurion (°) 
is used to enforce the AT = ^ rule for isospin violation, observed in the free A decay |18| . For 
the rho meson, we have 



n/W 

n ANp 



G F /4 tpN 



A p Yls + B v p Y + B T p 



a^d„ 



2M 



n 



NNp 



9nn p 7 + 9nn p ■ 



a 



2M 



(36) 



(37) 



The corresponding Hamiltonians for the r\ and uo are completely analogous to those of the 7r and 
p, respectively, if one takes into consideration their isoscalar nature. 

The weak couplings of the heavier mesons are theoretically inferred from those of the pion 
through unitary-symmetry arguments and other relationships. The strong ones are again taken 
from OME models for the nuclear force. We will follow the parametrization adopted in Ref. |18j . 
where further details can be found on this matter. For definiteness, the numerical values are 
reproduced in Tabled 



2.3.1 One pion exchange 

The local nonrelativistic one-pion-exchange transition potential is given in momentum space by 

= -{ 1 + %Y Gf a w - ■ - C 4 * + m - ' «) • (38) 

where 

M = EpM± M , (39 ) 

3M + M A v ; 

Comparing Eg . (J38|) with the result that would have been obtained under approximation (|31|) . 
namely, [H Eq.(24)} 

<40> 
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Table 1: Coupling constants, masses (/Xj) and cutoff parameters (Aj) for the nonstrange mesons 
The weak couplings are in units of Gf [J%- Adapted from Ref. 18 . 



Meson 






Coupling Constants 




Hi 


A. 


i 






Weak 


Strong 


[MeV] 


[GeV] 






PV 


PC 








7T 




= 1.05 


B n = -7.15 


Qnnit = 13.3 


140.0 


1.30 


V 




= 1.80 


B v = -14.3 


QNNr, = 6.40 


548.6 


1.30 


p 




= 1.09 


B v p = -3.50 
Bj = -6.11 


9nn p = 3-16 
9nn p = 13.3 


775.0 


1.40 


to 


A u 


= -1.33 


= -3.69 
Bl = -8.04 


9nnw = 10-5 
9nnlu = 3-22 


783.4 


1.50 



it is possible to estimate the relative size of the effects of the more accurate treatment of the 
kinematics, adopted here, from the following correction factors: 



(1 + M /M q 



0.998, M/M = 0.996 



(41) 



and 



fa/Hn = 0.752. (42) 

If each of these values were equal to unity, there would be no effect at all. Apparently, the 
situation is not much different from this, except in the last case, which will have a noticeable 
effect since it increases by ~ 35% the range of the corresponding contribution to the transition 
potential. 

When changing to the coordinate representation through Eqs. (|24j) or ((25)) . the shape func- 
tions 



fc(r, M) 
fv(r,p) 
fs{r,fi) 
fr(r, fi) 



a iqr 



d 3 q 

(2vr) 3 " 



1 



-fj/r 



q 2 + /x 2 



Airr 



1 

3 

1 2 

5" 



M 2 fc(r,fj) - 5(r) 



1 + — + 



fc(r,fi) 



(43) 



fj,r (pr) 2 

naturally arise, accordingly as the numerators in the Fourier transforms are, respectively, a 
constant, a vector, a scalar or a tensor built, at most quadratically, from q. In terms of these, 
we get, for the potential (|38|) in coordinate space, 



K (0) (> 



M n 



1 + wJ ^ F 



2 9NNn 



2M 



%A % f v {r,fl n )(T2 ■ r 
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B B 

-r^j. fs(r,fin)cri ■ <x 2 + -Z- f T (r, p, n ) S 12 (r) 



(44) 



where r = r/r and Si 2 (r) = 3(o"i • r )(cr 2 • r ) — a\ ■ a 2 . Under approximation (J31|) . this would 
reduce to 



V<°>(| 



n 2 9NNtt 

Gf ^ W T1 " T2 



iA„ f v (r, (t 2 r 



B B 



(45) 



The first-order nonlocality coefficient in momentum space, appearing in Eq. (j21j) . is given, 
for the pion, by 



, . Mi | ~ 2 9NNn (Bit \ &2 ■ q 



M„ 



2M 



V2M 



Q 2 + /4 



where 



11 1 

M ~ M ~ Ma ' 

Changing to the coordinate representation according to Eq. (1251) . we get 
\ADf \ f-i i M °\ n 2 9NNn B v _ 



2M 2M 



(46) 
(47) 

(48) 



and introducing this into Eq. (J28I) yields for the first-order nonlocality correction 

-i 



1 + 



Mo 
M„ 



^ 2 <7JVMr -Stt 
Cxi? ^ 



2M 2M 



— Tl • T 2 X 



\ 3M~ + M t^'^ " 1 ' 0-2 + Mr,p.„) Sn(r)] 



fv(r,fin) (<T2 ■ r) (<Ti • V)| . 



(49) 



The mass averaging approximation (J31J) would set 1/M to zero. Therefore, \q) = and 
there would be no first-order nonlocality correction for the pion under this approximation, i.e., 



^(i) 

VL } r =0. 



(50) 



2.3.2 One rho exchange 

The one-rho-exchange contribution to the local nonrelativistic transition potential in momentum 
space is 



1 + w) g fA t ^ t * [ K l- K U 2 



K P (°"l x Q) • (""2 x q) 
iKUax x a 2 ) ■ q + K b (tr x ■ q) 



q 2 + jj? D 



(51) 
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where, for notational convenience, we have introduced the coefficients 



K 



K 



Ki 



a p 9nn p 



tdV „V 

d p 9nnp 



(±S + (±.W l -(± 



+ 



\4M / \4M J 2 I M, 



/ Bp 9nn p Bp_ 9nn p 
\2M 2M 2M 2M 



1 + * 

M„ 



+ 



Bj 9 T NNp ( Ml \ 
2M 2M \AMM ) ' 



2M 2M 



V T 

9nn p 9nn p 



2M 



2M 



A, 



V T 

9nn p 9nn p 



2M 



2M 



1 + 



Mo 



A 9nn p (M \ 
p 2M \2MJ ' 



(52) 



The corresponding potential under approximation Vp°\q), can be obtained from 

Eq. (|51|) through the substitutions: 4 



' P y O 1 



K ] 
p 



Kl (j = l-5), M /M q 



ftp y ftp i 



(53) 



with 



K 



K 



Ki 



Kl 



D p 9NNf 
a p 9NNr 



\4M J 



' B v p+B y 

2M 



2 /J^ • 

+ \AM 
V I „T 



_ -°p ffiyTVp -°p fljVTVp 

2M 2M 2M 2M 



9nn p + S'atat, 



2M 



-4,, 



0. 



' V i T N 

2M 



(54) 



Let us now compare Eqs. (|52l) with the corresponding expressions under approximation (pTT 
namely, Eqs. (|54j) . Firstly, we notice that the two correction factors in Eq. (J41I) as well as 



fip/ tip = 0.992 

are very close to unity. Secondly, we also notice that the relative values of the remaining 
correction terms can be estimated from 



- - 

2 I M n 



{am) + \AMJ 



1.85 x 1(T 3 , 



4 This agrees with Eq. (34) of Ref. ^8], except for a wrong sign (A p — > —Ap) and an omitted term (oc q 2 ). 
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4MM ' 2M 

and are, therefore, considerably less than unity. We, thus, expect that, as far as the local 
contributions are concerned, only very small corrections will result from the more accurate 
treatment of the kinematics in the present case. 

Making use of Eq. (|24|) . we get, for the potential (|51|) in coordinate space, 



(l + G F £ n • r 2 [K] f c (r, jx p ) + 3K 2 p f s (r, fi p ) 

+ 2K 3 p f s (r, fi p ) ct\ • <x 2 K 3 p fr(r, fi p ) S 12 (f) 
+ fv(r, fip) K A p (<n x <x 2 ) + iK 5 p a 1 • r } . 



(55) 



Once more, the corresponding potential under approximation i\'Al\) can be obtained from the 
above equation by means of the substitutions (|53|) . 

For the rho meson, the coefficient of the linear term in Eq. (|21l) is given by 



where 



M n 



-i 



1 + WJ Gf^t.-tzx 



K p q- iKp <ti x q - iK p <x 2 x q 

- K 9 (ax x q) x er 2 + Kl° (<x 2 x q) x a x 



+ K)} <ti - iKf <Ti x er 2 



q 2 + fi 2 p 



(56) 



Kl 



Kl 



K 



10 



n P 9nn 



J_ ( A _ J_ i r 

p 16M \M M ' ' 



T 

v 9nn p 
p 



M 



+ 2M 9nn > 



1 

2M 



1 + 



Mo 
M r , 



2M 4MM 
M 



a p 9nN(. 



1 



2 j 1 
8M \M + M ] + 



2MM 
1 



+ 



Bp 9nn p 



M 



o 



8MM 



+ 5. 



2M 2M 4MM 
M 



T 

v 9nn p 



+ ^ SiVTVp 



1 



1 



M M 
1 

M 



M, 

T 

y ffiViVp 



+ 



2M 4MM 



M 2 



2M 2M 4MM 



2M 



1 



2M yAWp 4MM 2M 2M 2 



5 



M Bj g T NNp M 
p 2M 2MM 2M 2M M 



T 

v 9nn p 



1 + 



2 

11 

M \ 

M 



1 



'1 / M 
M \ 



M 



M„ 



1 + 



M„ 



+ 



M 
MM 



-°p 9 N N £ 



1 



4MM 
M 



+ 5 



T 

v 9nn p 



1 



+ ^I V 

2M yNNp 4MM 2M 2M 2M 



p 2M 2M 



1 + 



M q 

Mo 
M n 



B J 9nn p Mq 



1 + 



Mo 
M„ 



12 



Ap 9nn p 



1/21 
2 \M + M 



K 



12 



a 9nn p Mp 
p 2M M ' 



(57) 



with 



111 



(58) 



To get the first-order nonlocality correction Vp (r), we need first to change (|56|) to the coordi- 
nate representation, according to Eq. (I25j) . This gives 



V«(r) 



1 + TT G F /J I 'T 1 -T 2 X 



M, 



f fv(r, Ap) 



iK p b r + a-i x r + K° cr 2 x r 



iKp (<Ti X r) X <T2 + «^p° (<7"2 x r) x <Ti 



+ / c (r, Ap) o-i - ^p o-i x <t 



(59) 



Introducing (j59|) into Eq. (j28jl and noticing that — icr x r • V = <x • Z, where Z = — ir x V is the 
relative orbital angular momentum, we obtain, finally, 



1 + 



( 2M A 
1 3 Ma + M 



G F /4 Tl • T 2 X 



/ s (r, Ap) (3i^ 6 + 2(K 1 p ° - K 9 p ) <j x • <x 2 

10 7^9 



p 1 -v—p 
(K p lu -^ p y )/ T (r,Ap)S 12 (f) 



- /y(r, Ap) (i^f °"l x °"2 + i-K"" o-i ) • r 
fvijr, Ap) 



p 1 z 1 p 
Klr-V + Kla l -l + Kla 2 -l 



+ Kla 2 -r<r 1 -V-Kfa l -ra 2 -V 

+ / c (r, A P ) (^p 12 <tiX<t 2 + iKf a x ) ■ v| 



(60) 



Under approximation (|31|). the only surviving coefficients for the nonlocal potential would 



be 



Kl 



1 / 1 



1 / 1 



+ 



B P 9nn p iM ( m + m ) + ^ 9nn p [ 77 + 77 



2M 



1 1 

/p \ M + M 



K 



ii 



A V 

a p 9nn p 



1 j 1 

M + M 



, 9 N N P 



1 1 

+ 



2M \M M 



(61) 
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and Eq. (|6()j) would reduce to 

viV) = - 



Gf a4 t~i ■ r 2 



fv(r, Hp) 



(K 7 pt T 1 -l + K 8 p a 2 -l 



iK' L f c (r, n P ) a t ■ V + -iff 1 f v (r, p P ) <r x ■ r 



(62) 



2.3.3 Extension to 77 and w exchanges 

If one remembers to make the replacement T\'T% — > 1, the results obtained above for the ir and 
p mesons can be straightforwardly extended to the r/ and u>, respectively, and the corresponding 
expressions need not be reproduced here. Let us just mention that the ratios jl^f 'fj,^ = 0.985 and 
fiu/fJ-u = 0.992 are very close to unity and, consequently, as happened for the p, the effects of 
the reduction of the effective mass are much less important for these mesons than they are for 
the pion. 



2.4 Contributions of strange mesons 

For the strange mesons, the weak and strong vertices in Fig. n are interchanged with respect to 
those for the nonstrange ones, i.e., 



a = W, b = S (nonstrange mesons), 

a = S, b = W (strange mesons). 

For the kaon, the effective Hamiltonian for the strong coupling is |18l Eq.(28)\ 

Wank = ^9ank i>Nls ~4>A, 
while, for the weak one, it is [TBI Eq.(29)] 



n_,W 
n NNK 



iG F fi i/j n 



C P K V - 



+ 75 



IpN- 



For the K*, we have, for the strong coupling 18, Eq.(38)], 



'Hank* = 
and for the weak one ^Hl Eq.(39)], 

G F pl i) N {7 



v 



9ank* 7 + 9ank* 



2M 



n_/W 

n NNK* 



n pcy 



* } \t /0 N 



+ 



2M 

+ 7^75 



c k y(°\ (^ ) ) t +^(^ ) ) t f^ 



(63) 
(64) 



(65) 



(66) 



(67) 



We again follow the parametrization adopted in Ref. ^H], and collect the numerical values in 
Table 121 for convenience. 
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Table 2: Coupling constants, masses and cutoff parameters (Aj) for the strange mesons. 
The weak couplings are in units of Gf Adapted from Ref. 18 . 



Meson 




Coupling Constants 




Mi 


A* 


i 




Weak 


Strong 


[MeV] 


[GeV] 




PV 


PC 








K 


C? v = 0.76 
= 2.09 


C P K C = -18.9 
D PC = 6.63 


9ANK = -14.1 


495.8 


1.20 


K* 


CgY = -4.48 
DftY = 0.60 


C^ V = -3.61 
C£?' T = -17.9 
D™ ,V = -4.89 
Dk C *' T = 9.30 


9%K* = - 5 -47 
9ank* = -H-9 


892.4 


2.20 



These mesons are isodoublets, and in terms of their different charge states we can write 



1, 



to 



(^ ) ) t =(^ +) ) t - + ^(^ 0) ) t ( 1 -^ 



(68) 



for the kaon, and similar equations for the K* . As a result, when applying the Feynman rules 
to compute the transition potential as explained in Subsection 12.11 the isospurion will permit 
the introduction of isospin operators of the form 



1= - C(l+T 1 -T 2 )+D 



(69) 



for each of the different couplings in Eqs. (|65j) and (jfi7|) . Explicitly, they are 

1 
2 



I 



PV 
K 



Cf(l+T r T 2 ) + Df, 



I 



PC 
K 



1 



c£ c (l + n ■ r 2 ) + d 



PC 
K ) 



(70) 



for the kaon, and 



r PC,V 


1 

2 


C^,' (l + Tl-T 2 ) + %' 


T PC,T 
i K* 


1 

2 


C^,' (1 + Ti • T 2 ) + 


T PV 


1 

2 





(71) 



15 



for the K*. It then becomes apparent that each such coupling will give a contribution propor- 
tional to ^ C + D to the isoscalar potential and a similar one proportional to \ C to the isovector 
potential. 

2.4.1 One A' exchange 



The contribution to the local nonrelativistic transition potential in momentum space due to this 
meson is 

W c \ ai-q 

a2 q] ¥TW K - 



v (o) (n \ - (1 + M A 1 r,„ 2 9ANK (t pv 
Vk {q) -{ 1 + W q ) G ^W( /J( "2M 



(72) 



Comparing this with the result that would have been obtained under approximation (Jj2 
namely 5 , 



f>(°)/ \ n 2 9ANK ( t pv i k C \ °"1<? 



(73) 



and noticing that the two correction factors in Eq. ([41ft as well as 

fi K /fi K = 0.982 

are very close to unity, one can see that, for the kaon, only very small corrections will result in 
the local contributions from the more accurate treatment of the kinematics. The expression for 
the potential (|72[) in coordinate space is 



Vf(r) 



1 i M I n 2 9ANK 



2M 



jPC jPC 

f s (r, fi K ) °i-°2 + ^ Mr, jx K ) S 12 (r) 



(74) 



and, under approximation ([31 ft . this becomes 



vf (r) = G F £ 



2 9ANK 

2M 



jPC jPC 

+ -^fs(r,n K )<Ti'<T2 +^fT(r^ K )S 12 (r) 



(75) 



Starting with the first-order nonlocality coefficient in momentum space in Eq. ([21[K we have, 
for the kaon, 



1^(1)/ \ [-, , M o\ n 2 9ANK ( t pv i k° ^ 



q 2 + &' 



(76) 



which, in the coordinate representation, becomes 



V£>(r) 



1 + tt) 



M 
tPC 

2M 



2M 



fv(r,(lK) (<T2 ■ r)cri 



!k V fc(r, fix)*! 



(77) 



5 This differs from Eqs. (24) and (31) of Ref. |18) in the sign of and the interchange of <j\ and <Ji. 
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Introducing this into Eq. ()28j) yields, for the first-order nonlocal potential, 

-l 



V M q J 2M 

{ 9 Ma \T PC T pc 
} *-/s(r,/2*>i-*2 +4M^if)^(r) 



\ 3M\ + M [ 2M 

+ Hk V fv(r,p, K ) (en • r) 



2M 



tPC 
1 K 

2M 



Mn/fe-)fo -v) . (78) 



As already stated, the mass averaging approximation (J31|) would set 1/M, defined in Eq. (|47|) . 
to zero. Therefore, there would be no first-order nonlocality correction for the kaon under this 
approximation, i.e., 

V ( K(r) = 0. (79) 

2.4.2 One K* exchange 



For one--fT*-exchange, the local nonrelativistic transition potential in momentum space is 

-i 



Mr 



1 + wJ *~* F 



k\, - kh q 2 



- Kk* (<ti x q) ■ (er 2 x q) 

- i K K , (<ti x a 2 )-q + K%, (er 2 ■ q) 



1 



q 2 + //if. 



where we have introduced the isospin operators 



v jPcy 
Sank* a k* > 



R- 2 



V TJTI. 

Sank* l k* 



pcy 



+ 



1 / 1 



Um) + Um) ' 2 VAf, 



+ 



/ T r^C*,V y T P< 

I Sank* 1 k* _|_ Sank* 1 k 



SA K » 



n K* 



k 5 

JA K , 



\ 2M 2M 2M 2M 

I sInk* I%' T ( Ml \ 
2M 2M \4MM J ' 



1 + 



Mo 

Mn 



T PC,v t pc,t 

K* _|_ i K* 



I 



2M 

PV 

K* 



2M 



1 + 



V T 

Sank* , Sank* 
2M 2M 



1 + 



Mo 



v t 
Sank* , Sank* 

2M 2M 



K * 2M V2M/ 



(80) 



(81) 



The corresponding potential under approximation (|31|). V K ^)(q), can be obtained from 
Eq. (|8U|) through the substitutions: 



V, 



(o) 



A"* 



(0) 



K* ! 



(j = 1-5) , Mo/M q , /ijf. /xx* 



(82) 
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with 



K 2 



K 3 



K 4 



v jPcy 
9ank* 1 k* 

v jPcy 
9ank* l k* 



V4M J 



4M J 



+ 



T T^L 

9ank* 1 k* 



2M 2M 



+ 



9ank* a k* 



2M 2M 



' t PC ' V _i_ r PC ' T ' 
1 K* +1 K* 

2M 



jPV I 9ank* +_9ank* 
K * \ 2M 

0. 



9ank* + 9ank* 



2M 



(83) 



By an analysis very similar to the one performed for the p meson and noticing that 

Ak* / vk* = 0.994 

is also very close to unity, one concludes that, again in the present case, only very small correc- 
tions will result in the local contributions from the more accurate treatment of the kinematics. 
For completeness, we give below the expression for the potential (|8U|) in coordinate space: 

4°*V) = { l + Tf) {Khfc(r,PK*)+3Khfs(r,flK*) 



+ 2K\, f s (r, fi K *)<rx ■ CT2 ~ K 3 K * fr(r, fa.) S 12 (r) 

+ fv{r, m*) 



Kk* {(T\ x er 2 ) + iK b K * a 



K* 



(84) 



Once more, the corresponding potential under approximation (|M1I) can be obtained from Eq. (|84j) 
by means of the substitutions (|82|) . 

The first-order nonlocality coefficient in momentum space, appearing in Eq. (|21|) . for this 
meson, is 



- Kx* (a i X q) X <r 2 + Kk* (<x 2 x q) x <j 1 



Kk* <x 2 + iKx* <ri x cr 2 



q 2 + P? K * ' 



(85) 



where 



K 



K* 



K 



K* 



v t pc,v 

9ANK* i K* 



1 



1 



16M \M M 



v 



Tv _ 77 + ^AiVif 



/£F' T M 



2M 4MM 



+ 2M 



1 

2M 



1 + 



Mo 



Mr 







V T PC,V 

9ank* l k* 



1 



2 1 

8M VM + m. ' + 



2MM 
1 



+ 



9ANK* i K* 



PC,T 



Mn 



2M 2M 4MM 



8MM 



-4- n y ±K 

+ 9ank* 



M 



2M 4MM 



2iW ^* O i Jl/f nV ; I 71 /T I 71/T 1 



2 \M M 



PC,T 



Ml 



2M 2M 4MM 
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K 



K * 



v jPcy 1 ( 1 , 1 
9knk**k* w{~M + ti 



V 

9ank 



jPC,T 1 



2 j 1 
2M 2 VM + M~ 



1 + 



Mo 
M„ 



, ffAjVX* rPC,V 

+ 2M 



M< 



o 



4MM 



+ 



T T PC,T j. j, 

9knk* i k* m o 
2M 2M 2 



1 

M 



1 + 



Mo 
M„ 



+ 



M 
MM 



k 9 



k 10 

K* 



V A K 
9ANK 



jPC,T ; t>~^- 

I Mo_ 9ank* L K* 



T T PCT M 



2M 2MM 



— o / ^ Mq 
2M 2M M \ M q 



9ank* 1 k* 



pcy 



1 



4MM 



4- n y Ir 

+ Sanx 



PC,T 



1 



M 



2M 



A' 



ii 

K* 



4MM 
1/2 1 



- + 



^* J ^2lF + Sj + 2M 



2M 2M 

T T PC,T 

2M 2M 2M 
9'ank* t pv m o 

1 If* 



1 + 



Mo 
M„ 



1 + 



Mo 

M r , 



2M 



K 12 



V 

9ank 



tPV 1 , 9ank* t pv 

Mq 

K * 2M + 2M lK * 2M ' 



(86) 



with 1/M and 1/M as defined in Eqs. (|47j) and (|58|) . To get the first-order nonlocality correction 
vj()(r), we need first to change (|85j) to the coordinate representation. This gives 



V£l(r) 



1 + 



Mo 

M n 



-1 



{ fv(r, fix* 
\ r 



G F f4 x 



iK%* r + K 7 K * ai xr + k^* a 2 x r 



iK\* {cr\ x r) x <r 2 + iA^4 (er 2 x r) x er 



er 2 - o"! x er 2 



Introducing (|87|) into Eq. (J28|) . we obtain, finally, 



(87) 



-(^-KiO/r^/i^)^^) 
+ fv{r, m*) {k]t, a 1 X er 2 + cr 2 ) • f 



jV(r, £x* 



r ■ V + i^* a l ■ I + Af , c 2 • Z 



<ti • cr 2 r • V 



+ A^, <x 2 • r <ti • V - a 1 ■ r a 2 ■ V 



- /o(r, Ax*) (Kg. a 1 xa 2 + ik}, 1 , <x 2 J • V 
If one assumed that the averaged-mass approximation (|31|) could be made, several terms in 
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the nonlocal potential would disappear. The only remaining coefficients would be 

K 7 n v T PC ' V 1 ( 1 + 2 ^ 4- 9 ^ NK " T PC ' V ( 1 + 1 

" 9ank* 1 k * aM \ M + M J + ^* \Af M 

V 8 _ V T PC,V 1 f 1 , 2 \ , „V J K* ' T ( 1 1 



Kk * ~ Sank* *k* (jg + ^ 
and the first order nonlocality correction would reduce to 



2 



fv(r, t*K*) 



^ (k k * <T 1 ■ I + R K . <T 2 ■ I 



— 7 — 

+ iKft* fc(r, Hk*) cr 2 • V - -if^ 1 * jy (r, \i K * )cr 2 -r 



(90) 



It is interesting to point out that, whereas the first-order nolocal terms are systematically 
omitted in the literature on nonmesonic decay, they are routinely included in the closely related 
domain of strangeness-conserving, parity-violating nuclear forces. (See, for instance, Eq. (115) 
in Ref. We note, however, that the terms proportional to K 11 and K , for the vector 

mesons, have been recently discussed in the literature |21| . 

2.5 Finite size effects 

Before closing this section, let us mention a refinement that should always be added to the strict 
OME description we have been developing up to now, especially when large momentum transfers 
are involved, as is the case for nonmesonic hypernuclear decays. This is the effect of the finite 
size (FS) of the interacting baryons and mesons at each vertex. 

Taking a clue from the OME models for the NN force 38, 39 , the FS effects are phenomeno- 
logically implemented in momentum space by inserting, at each vertex in Fig. ^ a form factor, 
which we choose to be of the monopole type, 

A '~^ (91) 

where i refers to the meson involved and Aj are the cutoff parameters in Tables ^ an d El This 
corresponds in coordinate space to replacing, in the expressions for the transition potential 
discussed in Subsections 12.31 and 12.41 each of the shape functions (|43|) as follows: 

fN(r,fii) - f N (r,p,i) - f N (r,Ai) + ^-/^(r, AJ, (92) 

where N = C, V, S, T. When the kinematical effects are ignored, Eqs. (19 ip and (|92l) should be 
modified by making /2j — ► /ij, thus leading to agreement with Ref. [T8| . 

In what follows, it is to be understood that these FS effects are always included. 

3 Numerical results and discussion 

3.1 Decay rates 

We present here the numerical results for the different contributions to the nonmesonic weak 
decay rates of 1 ^C. We consider, separately, the neutron-induced (n) and the proton-induced 



20 



(p) contributions, as well as those coming from the parity-conserving (PC) and parity- violating 
(PV) transitions. All quantities are in units of the free A decay constant, Tq = 2.50 x 10 -6 eV. 
The main observables are the total nonmesonic decay constant T nm = T n + T p and the ratio 
T n /T p , whose experimental estimates are in the ranges 0.89 - 1.14 and 0.52 - 1.87, respectively, 
with large error bars |46j-|50j. Most, if not all, calculations in the context of OME models give 
reasonable results for T nm but fail completely for T n /Y p . However, our main objective here is 
not so much to try to reproduce the experimental values for these observables, but rather to 
assess the relative importance of the kinematical and nonlocality effects, usually ignored, in their 
theoretical prediction. For simplicity, we restrict the discussion of the nonlocality corrections to 
those of first-order. 

For the explicit evaluation of the transition rates, we follow the approach of Ref. 21\ . The ini- 
tial and final nuclear states in Eq. (^Q) are described in the extreme particle- hole model (EPHM) , 
taking as vacuum the simplest possible shell-model approximation for the ground state of 12 C, 
namely, lsi/2 an d IP3/2 orbitals completely filled with neutrons and protons. 6 The A single- 
particle state has quantum numbers j\ = lsi/2 and the nucleon inducing the transition occupies 
a 32 = 1^1/2 or I-P3/2 orbital. Therefore, [5J Eqs. (4-2,3)] 

\I) = KM) (jn)- 1 ; Jj) , \F) = l(jn)- 1 (j 2 N)~ l ;J F ) , (93) 

where «// = 1, j = lp^/2, N = p or n for proton- or neutron- induced transitions, respec- 
tively, and J F takes all the values allowed by angular momentum coupling and (when relevant) 
antisymmetrization. Then, changing the momentum variables in Eq. (^Q) to relative (p') and 
center-of-mass (P') momenta, making a multipole decomposition of the corresponding free waves 
and performing the angular integrations, one gets, for iV-induced transitions, [2"TI Eqs. (2.4,9)] 



16M 3 

N z 



vr 



£ / de'y/e>(A hN -e>) 



\^'P'L r X'S , J l T'M! r ,(jnr 1 U2Nr 1 ;J F -,J I \Hh^(3nr 1 -,Ji) , (94) 



I'L'X'S' 
J'T'M' 



where T' is the total isospin of the two emitted nucleons and the angular momentum couplings 
V + Li = A', A' + S' = J' and J' + J F = J 1 are carried out. One also has P' = 1\J Me', 



p = yM(Aj 2 jv — e') and Aj 2 n = Ma — M + Ej 1 a + £j 2 N, where the single-particle energies are 
taken from experiment, according to Table 3 of Ref. |15j . 

After some standard manipulations, Eq. (j94jl takes the form |21| Eqs. (2. 13), (4-4)] 

16M 3 



vr *-r" Jo v 



n 

Y, F J> N E \M{p'l'P , L'>!S'J , T'-j 1 K32N)W (95) 

J' I'L'X'S'T' 

which allows a nice separation between the nuclear structure aspects and those of the decay 
dynamics proper. The nuclear structure factor is, in second-quantized notation, 



1 Jf 



6 As shown in that reference, further sophistication of the nuclear structure description has little effect on the 
nonmesonic decay rates. 
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and its nonzero values, for the nuclear states in Eqs. (|93|) . are: F Q Sl/2n = F Q Sl/2P = 1/2, 

F ls 1/2 n = F ls 1/2P = F l P3/2 n = F l P3/2 n = ^ p l Pz/2P = ^ ^ p l Pz/2P = ^ ^ 

Table I in Ref. [131 ■) O n the other hand, the nuclear matrix element governing the decay is 

Mip'l'P'L'X'S'J'T' ■ j t A j 2 N) = 

( P 'l'P'L'\'S'J'T'M T \ V \jiA j 2 N J') , (97) 

where (• • -| V |- • •) is a direct matrix element and the factor in front takes care of antisymmetriza- 
tion. To compute the isospin part of this matrix element, one writes the baryon content of the ket 
as \AN) = \^m tA ^mt N ), where m tN takes the values m tp = 1/2 for protons and m tn = — 1/2 
for neutrons, while, in accordance with the isospurion stratagem, one treats the A as if it corre- 
sponded to mf A = —1/2. On the bra side, one sets Mt = mt A + m tN . To simplify the spatial 
integration, one resorts to a Moshinsky transformation |5J of the initial AiV system. To this 
end, the shell-model radial wave functions are approximated by those of a harmonic oscillator 
with a length parameter of b = 1.75 fm, which is an average between the values appropriate for a 
A and for a nucleon [EJ. Some useful expressions for the computation of these matrix elements 
are given in Appendix lAl 

As done in Ref. |21j and already stated above, the FS effects are taken into account as 
indicated in Subsection 12.51 Another important effect to include due to the relatively large 
momentum transfers involved in nonmesonic decays is that of short range correlations (SRC). 
The most satisfactory way to deal with the SRC between the A and the inducing nucleon in the 
initial state would be through a finite- nucleus G- matrix calculation |52) . However, as mentioned 
in Ref. JE]> this can be well simulated by means of the correlation function 

9 AN(r) = (l - e"^ 2 ) 2 + fir 2 e~ r '' ^ , (98) 

with a = 0.5 fm, f5 = 0.25 fm~ 2 and 7 = 1.28 fm. As for the SRC between the two emitted 
nucleons, one might want to perform a T-matrix calculation including final state interactions 
along the lines of Ref. 53 .'' A simpler, if less satisfactory, way is to again appeal to a correlation 
function, like 54 

9NN{r) = 1 -jo(q c r) , (99) 

where jo is a spherical Bessel function and q c = 3.93 fm -1 . For our purposes here, it is sufficient 
to follow Ref. [U and opt for these phenomenological correlation functions. Thus, in the 
calculation of the nuclear matrix elements in Eq. f|97l) we simply make the replacements 

\jiAj 2 NJ) -» g AN (r)\jiAj 2 N J) , 
(p'l'P'L'X'S'j'T'Mfrl -► (p'l'P'L'X'S'J'T'M^\ g NN {r). (100) 

Following the discussion in Subsections 12.31 and 12.41 we initially focus our attention on the 
reduction of the effective meson masses, jli in Eq. f|20j) . especially that of the pion, jl n in Eq. ()42l) . 
and show that indeed this is the main kinematical effect for the local potential, but not so for 
the nonlocal one. To this end we give, in Tables El and the corrections that should be added, 
according to several different calculations, to the standard OME results, i.e., those obtained 
when both the kinematical and the nonlocality effects are completely ignored. In Table 01 are 

7 In fact, there are claims that this is very important for a good description of the nonmesonic decay observables 
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Table 3: Corrections due to the kinematical effects on the nonmesonic decay rates of in 
several OME models, when only the local potential is included in the calculation. See text for 
detailed explanation. 



Model /Kinematics 


r PC 

1 n 


V PV 

1 n 


V PC 


r PV 
1 p 






7T 

averaged 

fllll 




0.0019 
n nm q 




0.0230 
D 0994 




0.0974 
n 1 nnn 




0.0601 
u.uoou 




0.1823 

D 1 S9Q 
u. loiy 




0.0031 
D DD94 


(vr,r/,K) 














averaged 
m -» fk 

full 




0.0024 
0.0023 
D 0(193 




0.0292 
0.0359 
D 0353 




0.0615 
0.0498 
D D5D8 




0.0653 
0.0707 
f) D694 




0.1583 
0.1586 
f) 1 577 




-0.0275 
-0.0150 

— n m 56 


vr + p 














averaged 

fJ-i -> fii 

full 




0.0008 
0.0008 
0.0008 




0.0214 
0.0209 
0.0206 




0.0712 
0.0711 
0.0729 




0.0658 
0.0685 
0.0652 




0.1593 
0.1614 
0.1597 




0.0055 
0.0047 
0.0046 
















averaged 
/j, n > [in 

IK -> fii 
full 




0.0048 
0.0058 
0.0058 




0.0245 
0.0323 
0.0306 




0.0710 
0.0605 
0.0619 




0.0803 
0.0953 
0.0878 




0.1807 
0.1940 
0.1862 




-0.0103 
-0.0032 
-0.0033 



the corrections corresponding to calculations that use only the local potential, and in Table |IJ 
those corresponding to calculations that include also the first-order nonlocality terms. In each 
table, the first column indicates which mesons have been included in the exchange process and 
how far the kinematical effects due to the lambda-nucleon mass difference have been taken into 
consideration. The entry "averaged" means that the mass-averaging approximation (J31|) has 
been made and, consequently, no kinematical effects have been included, while p^ — ► fj^ or 
Mi — * fti indicates that they partly have been, through these replacements made, respectively, 
for the pion alone or for all the mesons, in the expressions for the mass-averaged potentials V^ ^ 

and v\ ^ in Subsections I2.3l and l2. 41 (Excluding, of course, the factor Gf a4-) Finally, the entry 
"full" means that the kinematical effects have been fully taken into account, by making use of 
the complete expressions for V^ ^ and V^ 1 ^ when constructing the local transition potential and 
(for Table its first-order nonlocality correction. 

Examining the first block in Table one notices immediately that, when only the local 
potential is included in the calculation, the kinematical effects are well represented, for the pion, 
by the replacement fj, n — > p, n in the expression ([45 jl for the local transition potential obtained 
when they are completely negleted. Thus, the further modifications caused by these effects in the 
local potential, which lead to the "full" expression (|44|) . are of less importance. Comparing the 
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Table 4: First-order nonlocality corrections for the nonmesonic decay rates of in several 
OME models, and for different treatments of the kinematical effects. See text for detailed 
explanation. 



Model /Kinematics 


r PC 

1 n 


V PV 

L n 


V PC 


V PV 


r 

-i nm 




7T 

averaged 

fllll 





U.UUOi 





n 









n 
u 





U.UIUl 





— n nnfi9 

U.UUUi 


(vr,r/,K) 














averaged 
m -» fk 

full 






n 0037 






—0 0332 






031 2 






—0 0335 






— 031 8 






—0 0401 


vr + p 














averaged 

fJ-i -> fii 

full 


0.0001 
0.0001 
0.0001 
0.0023 


0.0016 
0.0017 
0.0018 
0.0014 


0.0004 
0.0003 
0.0003 
0.0386 


-0.0102 
-0.0110 
-0.0114 
-0.0088 


-0.0080 
-0.0088 
-0.0091 
0.0334 


0.0034 
0.0032 
0.0033 
-0.0001 
















averaged 
/j, n > [in 

IK -> fii 
full 


0.0013 
0.0014 
0.0015 
0.0133 


-0.0471 
-0.0500 
-0.0522 
-0.0901 


-0.0008 
-0.0011 
-0.0011 
0.0404 


-0.0970 
-0.1027 
-0.1071 
-0.1425 


-0.1435 
-0.1524 
-0.1589 
-0.1790 


-0.0238 
-0.0227 
-0.0229 
-0.0515 



last two lines in the remaining blocks of this table, one concludes that the analogous statement 
holds also for the other mesons. Finally, knowing this and comparing the second and third lines 
in these same blocks, one can see that the main local kinematical correction is that affecting the 
pion exchange. All these conclusions are in agreement with the discussion in Subsections I2.3.H 
I2T21 12331 HH and I2T21 

Going now to Table 01 one sees that the situation is quite different as regards the influence of 
the kinematical effects on the nonlocal potential. In fact, the first two blocks show that, in OME 
models containing only pseudoscalar mesons, the first-order nonlocality corrections vanish unless 
one takes the kinematical effects fully into account. This is just a restatement of Eqs. (|5()l) . ()79j) 
and the analogous result for the r\ meson. Similarly, examination of the last two blocks shows 
that, in OME models containing vector mesons, if one does not take the kinematical effects 
fully into account the first-order nonlocality corrections generally turn out very different from 
their actual values. The mere replacement pi — > jli does not work well in this case. This is so 
because, as can be seen in Subsections 12 . 3 . 2*1 and |2~4~21 several nonlocal terms appear as a direct 
consequence of the kinematical effects, rather than simply being modified by them. Therefore, to 
be consistent, one should take the kinematical effects due to the lambda-nucleon mass difference 
fully into account when dealing with the nonlocality corrections. 
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Table 5: Analysis of the different contributions to the nonmesonic decay rates of in sev- 
eral OME models. All corrections are computed with "full" kinematics. See text for detailed 
explanation. 



A/lnnpl / ( !fYn"i"T*"i rvn'hrvnci 

1V1UU.C1/ \_y vjllbl 1 U U U1U11D 


r PC 

1 n 


1 n 


l p 


l p 


r 

- 1 nm 


r iv 


7T 

Uncorrected value 
Local kinem. corr. 
1 st -order nonloc. corr. 

( r\ f vrp p~fr p n un lit p 


0.0063 
0.0019 
0.0032 


0.1162 
0.0224 


U. lOOU 


0.6019 
0.1000 
0.0729 
D 7748 


0.2887 
0.0586 


U.O^ / o 


1.0131 
0.1829 
0.0761 
1 9791 


0.1375 
0.0024 
-0.0062 
n i ^7 

U. -LOO / 
















Uncorrected value 
Local kinem. corr. 
1 st -order nonloc. corr. 

( .nvvpptprl unlit p 


0.0056 
0.0023 
0.0037 

n m 1 6 


0.2345 
0.0353 
-0.0332 
9366 


0.2124 
0.0508 
0.0312 
0.2944 


0.3813 
0.0694 
-0.0335 
D 41 72 


0.8339 
0.1577 
-0.0318 
9598 


0.4045 
-0.0156 
-0.0401 

f) 3488 


vr + p 














Uncorrected value 
Local kinem. corr. 
1 st -order nonloc. corr. 
Corrected value 


0.0056 
0.0008 
0.0023 
0.0087 


0.1004 
0.0206 
0.0014 
0.1224 


0.5140 
0.0729 
0.0386 
0.6255 


0.3608 
0.0652 
-0.0088 
0.4172 


0.9807 
0.1597 
0.0334 
1.1738 


0.1212 
0.0046 
-0.0001 
0.1257 


(7r,r l ,K) + (p,u J ,K*) 














Uncorrected value 
Local kinem. corr. 
1 st -order nonloc. corr. 
Corrected value 


0.0241 
0.0058 
0.0133 
0.0432 


0.2218 
0.0306 
-0.0901 
0.1623 


0.2961 
0.0619 
0.0404 
0.3984 


0.6238 
0.0878 
-0.1425 
0.5691 


1.1657 
0.1862 
-0.1790 
1.1729 


0.2672 
-0.0033 
-0.0515 

0.2124 



To better visualize our findings, we exhibit in Tableman analysis of the different contribu- 
tions to the nonmesonic decay rates of in the four OME models we have been considering. 
In the first line of each block, we give the values that would be obtained for the transition 
rates in the standard OME approach, i.e., when neither the kinematical, nor the nonlocality 
corrections are included. On the second line, are the corrections to these values arising from the 
kinematical effects related to the lambda-nucleon mass difference, but still restricted to the local 
contributions only. On the third line, we have the first-order nonlocality corrections, and on the 
last one, the values of the decay rates including the two corrections. As required by consistency, 
according to our previous discussion, both corrections are computed with the kinematical effects 
fully taken into account. 

Examining this table, one notices that the kinematical and the nonlocality corrections are 
typically of comparable sizes. Furthermore, for the partial and total decay rates, the former 
ones are always positive, while the latter are sometimes negative. Consequently, the two cor- 
rections should be included simultaneously, or not at all. Another point to remark is that the 
modifications in the uncorrected values of these decay rates when going from one OME model to 
another are of the same general magnitude as these corrections within each model. Therefore, 
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it might be questionable to consider other mesons besides the pion without, at the same time, 
including the kinematical and nonlocality corrections. 

The influence of the two effects together in the several partial decay rates varies around 
~ 80% for T PC and ~ 20% for the other ones, depending on the OME model. The net effect on 
the main decay obervables is smaller: it is ~ 15% for T nm and ~ 10% for the ratio T n /T p , again 
depending on the OME model considered. As one can see, these corrections are of no help to 
solve the discrepancy between the theoretical prediction and the experimental determinations 
for the latter quantity. They are too small for that, and usually go in the wrong direction. 
As a final observation, notice that the combined correction affects very differently the parity- 
conserving and the parity- violating transitions, especially when strange mesons are involved. For 
instance, in the model including the full pseudoscalar and vector meson octets (fourth block), 
the uncorrected value for Y PC /T PV is 0.379, while the corrected one is 0.604. 

One might wonder how important the corrections are for each meson-exchange taken in 
isolation. To answer this question, we show in Table El the different contributions to the partial 
and total decay rates, as well as to the n/p ratio, coming from each meson acting alone. Of 
course, in actual fact the contributions of the several mesons interfere with each other, so that 
the numbers shown in this table do not have any direct physical meaning, but they serve as an 
indication of the importance of the two corrections for each meson. 

The local kinematical corrections affect the pion more than any other meson, as expected 
due to its low mass. There, the effect on the decay rates is of the order of 20 - 30%, while for 
the other mesons it does not exceed ~ 10%. For the n/p ratio, the effect is always negligible. 

The first-order nonlocality corrections vary a lot, both in sign, and in magnitude. Depending 
on the meson and on the partial rate considered, the corrections can be as low as a few per cent, 
but in many cases reach the 50% level. The effect on the omega-exchange is exceptionally large, 
specially for the parity-violating transitions, where, in fact, more than half of the contribution 
comes from the nonlocal terms. However, for this meson the parity-conserving transitions are 
far more important and, besides this, the net nonlocal PC contribution is opposite to the PV 
one, such that the final effect of the nonlocal corrections on the total decay rate is of less 
than 5%. This is again an indication that the effects we are considering here are not exactly 
small, although, due to a series of compensations, they end up by not affecting the usual decay 
observables too much. 

3.2 Asymmetry parameter 

The only remaining nonmesonic decay observable, beyond T nm and r n /r p , that has been mea- 
sured to date is the asymmetry parameter, a a, which depends on the interference between the 
amplitudes for PC and PV proton-induced transitions to final states with different isospins. This 
parameter is a characteristic of the nonmesonic decay of a polarized A in the nuclear medium, 
having been defined so as to subdue the influence of the particular hypernucleus considered [T5] . 
It is experimentally extracted from measurements of the asymmetry, A y , in the angular distribu- 
tion of the emitted protons in the nonmesonic decay of polarized hypernuclei |55| I56j . There are 
large discrepancies, both experimentally and theoretically, in the determination of a\ [3J Sec. 7], 
specially after the newest experimental results for the decay of ^He ,56 a , which give a positive 
value for this observable, differently from all previous measurements. In strong disagreement 
with that, all calculations so far find a negative value for a\ [3], which makes the investigation 
of the kinematical plus nonlocality corrections on this observable particularly relevant. 

For the case of ^He, the following expression can be used to compute the asymmetry param- 
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Table 6: Analysis of the different contributions to the nonmesonic decay rates of A C coming 
from each meson acting alone. The interference among the different mesons is ignored. All 
corrections are computed with "full" kinematics. See text for detailed explanation. 



Meson / Contributions 


V PG 
1 n 


V PV 

1 n 






r 

L nm 




7T 

Uncorrected value 
Local kinem. corr. 
1 st -order nonloc. corr. 
Corrected value 


0.0063 
0.0019 
0.0032 
0.0114 


0.1162 
0.0224 


0.1386 


0.6019 
0.1000 
0.0729 
0.7748 


0.2887 
0.0586 


0.3473 


1.0131 
0.1829 
0.0761 
1.2721 


0.1375 
0.0024 
-0.0062 
0.1337 


V 

Uncorrected value 
Local kinem. corr. 
1 st -order nonloc. corr. 
Corrected value 


0.0020 


0.0003 
0.0023 


0.0031 
0.0002 


0.0033 


0.0047 
0.0002 
0.0005 
0.0054 


0.0024 
0.0001 


0.0025 


0.0122 
0.0006 
0.0007 
0.0135 


0.7200 
-0.0047 
-0.0130 

0.7023 


K 

Uncorrected value 
Local kinem. corr. 
1 st -order nonloc. corr. 
Corrected value 


0.0058 
0.0002 
-0.0002 
0.0058 


0.0453 
0.0041 
-0.0233 
0.0261 


0.0780 
0.0061 
0.0073 
0.0914 


0.0277 
0.0025 
-0.0139 
0.0163 


0.1569 
0.0129 
-0.0303 
0.1395 


0.4841 
0.0009 
-0.1886 
0.2964 


P 

Uncorrected value 
Local kinem. corr. 
1 st -order nonloc. corr. 
Corrected value 


0.0016 
0.0001 
-0.0005 
0.0012 


0.0017 
0.0002 
0.0007 
0.0026 


0.1274 
0.0070 
-0.0069 
0.1275 


0.0082 
0.0003 
-0.0013 
0.0072 


0.1390 
0.0075 
-0.0080 
0.1385 


0.0246 
0.0004 
0.0029 
0.0279 


Uncorrected value 
Local kinem. corr. 
1 st -order nonloc. corr. 
Corrected value 


0.0073 
0.0007 
0.0016 
0.0096 


0.0020 
0.0001 
0.0037 
0.0058 


0.0732 
0.0051 
-0.0030 
0.0753 


0.0019 
0.0001 
0.0011 
0.0031 


0.0843 
0.0062 
0.0034 
0.0939 


0.1232 
0.0028 
0.0708 
0.1968 


K* 

Uncorrected value 
Local kinem. corr. 
1 st -order nonloc. corr. 
Corrected value 


0.0062 
0.0001 
-0.0016 
0.0047 


0.0145 
0.0011 
0.0086 
0.0242 


0.0614 
0.0021 
-0.0025 
0.0610 


0.0189 
0.0007 
-0.0004 
0.0192 


0.1010 
0.0040 
0.0041 
0.1091 


0.2579 
0.0060 
0.0960 
0.3599 



eter [3]: 



a A = Ay&He) 



23? 



V3ae* -b(c* -V2d*) + V3f (\/2c* + d* 



|a| 2 + |6| 2 + 3 



|c| 2 + M 2 + |e| 2 + |/| 2 ) 



(101) 



27 



where 



a = (np/SolVIAp^So) 

= M(pl = PL = A = 5 = J=0T=1 M T = ;Ap, (ls 1/2 ) 2 .7=0), 

6 = ifap^PolVIAj^So) 

= iM(pl = l PL = A = l 5 = 1 J = 0T = 1 M T = ; Ap, (ls 1/2 ) 2 J = 0) , 

c = (np, 3 Si|V1Ap, 3 Si) 

= -M(pl = PL = A = 5=1 J=l T = M T = ;Ap, (ls 1/2 ) 2 J = l), 

d = -(rip/T^ |f|Ap, 3 Si) 

= M(pl = 2 PL = A = 2 5 = 1 J=l T=0 M T = ;Ap, (ls 1/2 ) 2 J = 1) , 

e = i(np, 1 Pi|V'|Ap J 3 Si) 

= -iM(p/ = l PL = A = l 5 = J=l T=0 M T = ;Ap, (ls 1/2 ) 2 J = l), 

/ = -i(np, 3 Pi\V\Ap, 3 Si) 

= -iM(pl = l PL = X = l S = IJ = 1T=1 Af T = 0;Ap, (ls 1/2 ) 2 J = 1) . (102) 

The extra factors in the transition amplitudes are due to differences in phase conventions, as 
explained in Appendix El and we have rewritten them in terms of the nuclear matrix elements 
defined in Eq. (1571) . 

It is important to realize that the formula (|101|) is not of general validity, and is, in fact, 
merely an approximation taken over from the result valid for a free space process |Uj and adapted 
somehow to the hypernuclear decay situation. In particular, there is no definitive prescription 
on how to devide the energy liberated in the decay, A_p, between the relative and CM motions of 
the emitted nucleons. It seems that, based on the expectation that the final result be insensitive 
to this point, some authors merely take P = in Eqs. (jl()2j) . while others integrate over phase 
space, both the numerator, and the denominator in Eq. (|101|) . We have checked that the two 
prescriptions indeed give similar results for ^He, and opted to tabulate only those corresponding 
to the first one. A more rigorous calculation of a\ is planned for the near future. 

We give, in Table [7[ the results we obtained for the decay rates and asymmetry parameter of 
A He. The calculation goes along similar lines to those explained in Subsection lH.il except for the 
obvious changes. The (hyper)nuclear model-space is restricted to the lsi/ 2 orbital; the relevant 

nuclear structure factors, Eq. fU), are |2H] F 1Sl/2n = F^ 1/2P = 1/2, F^ 172 " = F 1 1Sl/2P = 3/2; 
we used 1.62 fm for the oscillator length parameter; and took = 153.83 MeV. It is clear from 
this table that the comments made above about the decay rates of remain qualitatively valid 
in the present case. As to the asymmetry parameter, the effect of the two corrections reaches 
~ 18% in the ir + p model, but is of only ~ 5% in the complete model. On the average, it varies 
around ~ 10%. It is interesting to observe that the two corrections always go in the direction of 
making a\ even more negative, thus confirming the sign of the existing theoretical predictions 
for this observable in contraposition to its most recent experimental determination. 

4 Summary and conclusions 

We have proposed an approach that naturally establishes a hierarchy for the different levels of 
approximation in the extraction of the nonrelativistic transition potential in OME models in 
general. The central result is Eq. (|21j) . The first term corresponds to the local approximation, 
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Table 7: Analysis of the different contributions to the nonmesonic decay rates and asymmetry 
parameter of ^He in several OME models. All corrections are computed with "full" kinematics. 
See text for detailed explanation. 



Model /Contributions 


71 


71 


r£ 6 " 

p 


p 


1 


OA 


Uncorrected value 


0.0004 


0.0739 


0.3889 


0.1479 


1.7553 


-0.4351 


Local kinem. corr. 


0.0005 


0.0130 


0.0599 


0.0258 


-0.0295 


-0.0084 


1 st -order nonloc. corr. 


0.0011 





0.0461 





0.1810 


-0.0021 


Corrected value 


0.0020 


0.0869 


0.4949 


0.1737 


1.9068 


-0.4456 


{7T,rj,K) 














Uncorrected value 


0.0013 


0.1734 


0.1261 


0.2077 


0.3342 


-0.5852 


Local kinem. corr. 


0.0008 


0.0234 


0.0280 


0.0327 


0.0229 


-0.0106 


1 st -order nonloc. corr. 


0.0013 


-0.0285 


0.0176 


-0.0215 


0.0952 


-0.0406 


Corrected value 


0.0034 


0.1683 


0.1717 


0.2189 


0.4523 


-0.6364 


vr + p 














Uncorrected value 


0.0001 


0.0648 


0.3135 


0.1939 


1.2118 


-0.2665 


Local kinem. corr. 


0.0001 


0.0121 


0.0433 


0.0299 


-0.0251 


-0.0217 


1 st -order nonloc. corr. 


0.0004 


0.0016 


0.0245 


-0.0056 


0.1003 


-0.0273 


Corrected value 


0.0006 


0.0785 


0.3813 


0.2182 


1.2870 


-0.3155 
















Uncorrected value 


0.0112 


0.1672 


0.1778 


0.3642 


0.3558 


-0.5131 


Local kinem. corr. 


0.0026 


0.0206 


0.0344 


0.0442 


0.0233 


-0.0090 


1 st -order nonloc. corr. 


0.0054 


-0.0715 


0.0237 


-0.0897 


0.2073 


-0.0167 


Corrected value 


0.0192 


0.1163 


0.2359 


0.3187 


0.5864 


-0.5388 



usually adopted in the literature on nonmesonic decay. The second one, to the first-order 
nonlocality correction, which we have included in our calculations here. And the last one, to the 
second-order nonlocality correction, which we have neglected. We have also given a detailed and 
general account on how to deal accurately with the kinematical effects that result when one has 
different baryon masses on the four legs in the OME Feynman amplitude in Fig. ^ All this was 
particularized to A hypernuclear nonmesonic decay and detailed expressions for all contributions 
to the transition potential coming from the exchange of the complete pseudoscalar and vector 
meson octets were given. 

Using this formalism, we have investigated the relative importance of two effects sistemati- 
cally ignored in OME models for the nonmesonic weak decay of A-hypernuclei. First, that of an 
accurate treatment of the kinematics, i.e., of taking into account the difference in mass between 
the hyperon and the nucleon, when determining the OME transition potential. Secondly, we 
considered the influence of the first-order nonlocality-correction terms. Surprisingly, in view of 
the nonnegligible value of the mass-asymmetry ratio in Eq. (|32jl , we came to the conclusion that 
the kinematical effect on the local potential is small, except for the reduction of the effective 
mass of the pion, Eq. (|42|) . which implies an increase of ~ 35% in the range of the correspond- 



29 



ing transition potential. However its indirect influence is important, since it activates several 
nonlocal terms in the transition potential. 

Our conclusion is that the influence of the two effects together on the partial decay rates is 
sizeable, the full amount depending on which mesons are included. It can be very large for T PC , 
often exceeding 100%, while it typically stays in the 20-30% range for the other partial rates. 
The effects are somewhat washed out in the main decay observables, averaging to ~ 15% for the 
total nonmesonic decay rate, T nm = T n + T p , and to ~ 10% for, both the neutron- to proton- 
induced ratio, r n /r p , and the asymmetry parameter, a\. In particular, they do not in any way 
account for the well known discrepancy between the standard OME predictions for the latter 
two observables and their measured values [3]. To summarize, although the kinematical and 
nonlocality effects can be very important for particular transitions, they end up by not affecting 
the main decay observables too much. One can partly understand this from the following two 
facts: 

1. The most affected transitions are by far the neutron- induced, parity-conserving ones. How- 
ever those contribute very little to T nm and to T n /T p , and not at all to a\, which is a 
property of proton-induced transitions only. 

2. The largest relative effect of the nonlocality corrections is on the PV transitions coming 
from the one-omega-exchange process. However, for this meson, the PC transitions are 
much more important and they are affected in the opposite direction. As a result the 
two effects are largely neutralized in any of the main decay observables, which depend on, 
either the sum of the intensities, or the product of the amplitudes for these two types of 
transitions. 

We have also shown that the parameter T PC /F PV is strongly affected in OME models that 
include strange mesons. The relative correction is of ~ 30% in the (ir + r] + K) model and of 
~ 60% in the complete model. However, none of the three nonmesonic decay observables that 
have been measured up to now is sensitive to this ratio. 

Let us finalize by saying that, although the effects we have studied here can be considered 
small in view of the imprecision of the available measurements and the degree of uncertainty 
presently existing in the parameters of OME models for nonmesonic decay (particularly coupling 
constants), they are not altogether negligible. In fact, in many cases they appear to influence 
the theoretical predictions by roughly as much as the inclusion of other mesons beyond the pion 
in the exchange process. It seems, therefore, that they have a part to play in OME models 
for nonmesonic hypernuclear decays, specially if one takes into consideration that more detailed 
and accurate experimental data on this issue is forthcoming. 
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Appendices 



A Calculation of nuclear matrix elements 

We collect in this appendix some useful expressions for the calculation of the several matrix 
elements contributing to Eq. (|97j). As mentioned below that equation, the first step is to factor 
out the baryon content of the initial state, 

b'lA j 2 N J) = \j t h J) I AN) = \nihh n 2 l 2 j 2 J) \AN) , (Al) 

and perform a Moshinsky transformation [HJ f° r its space-spin part, 

h \ ji 



\nihji n 2 l 2 j 2 J) = y/(2j! + 1)(2 J2 + 1) ]T ^(2A + 1)(25 + 1) I h \ j 2 

AS [ A S J 

x \nlNL\SJ) {nlNL\\mhn 2 l 2 \) . (A2) 

nlNL 

A look at Subsections 12. HI and 12.41 shows that the operators involved are always of the general 
form v(r) fi(tri, cr 2 , r, V). All matrix elements are diagonal in the quantum numbers L and J, 
and the needed results are listed below. Some of them have already been given in Ref. [2^ but 
are repeated here for completeness. 8 

(p'l' P'L X'S'J\ v(r) \nl NL AS J) = 

5vi6 x , x 8 s >s (P'L\ NL) (p'l'\ v(r) \nl) , (A3) 

(p'l 1 P'L X'S'J\ v(r)ai ■ a 2 \nl NL XSJ) = 

SviSyxSs's [2S(S + 1) - 3] (P'L\ NL) (p'l'\ v(r) \nl) , (A4) 

(p'l' P'L X'S'J\ v{r)S 12 (r) \nl NL XSJ) = 

(-) L+J+l 8 sl 5 S 'i\l 120(21' + 1)(2Z + 1)(2A' + 1)(2A + 1) 



X 1' I r) W lAWNDWM^nl), (A5) 



(-f'+ s (p'l' P'L X'S'J\ v(r)a 1 ■ r \nl NL XSJ) 
(p'l' P'L X'S'J\ v(r)a 2 ■ r \nl NL XSJ) 



(-) L+J+l Jq(2V + l)(2l + 1)(2X' + 1)(2A + l)(2S' + 1)(2S + 1) 




x (P'L\ NL) (p'l'\ v(r) \nl) , (A6) 



8 There are some misprints in the formulas for the matrix elements in the published version of Ref.|21). These 
misprints, however, do not occur in its preprint version, axXiv:nucl-th/0011092 
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(p'l' P'L X'S'J\ iv(r)(<n x <r 2 ) ■ r \nl NL XSJ) = 

{-) L+J+S (Wsi + <Wso) V / 12(2/' + l)(2/ + l)(2A' + l)(2A + l) 
( V I I 







(V 1 Al 


j S> 


1 




\ I LI'] 


I A 


J 


?1 



(P'L\ NL) (p'l'\ v(r) \nl) , (A7) 



( p Y P'L A'5' J| u(r)<ri • V |nZ NL XSJ) = 
(p'l 1 P'L X'S'J\ v(r)a 2 ■ V \nl NL XSJ) = 



(-) L+J+1 y/6(2V + 1)(2Z + 1)(2A' + 1)(2A + 1)(2S> + 1)(25 + 1) 

( V 1 / \ J A' 1 X\ \ S' 1 5 \ / S' 1 S \ 
X 1, ) { I L I' j { \ \ \ ) \ X J X' j 

x (P'L\ NL) (p'l'\ v(r)d v {l', I; r) \nl) , (A8) 
(p'l' P'L X'S'J\ iv(r)(a 1 x <r 2 ) • V \nl NL XSJ) = 



iX' 1 Al 


j S' 


1 




\ / L V j 


I A 


J 


5) 



(-) L+J+5 (fe'ofei + 6s>i6so) \/l2(2l' + 1)(2/ + 1)(2A' + 1)(2A + 1) 

( V \ I 
x ^ 

x (P'L\ NL)(p'l'\v(r)d v (l',l;r)\nl) , (A9) 

(-) S ' +S (p'l' P'L X'S'J\ v(r)a- 1 ■ I \nl NL XSJ) = 
(p'l' P'L X'S'J\ v(r)a 2 ■ I \nl NL XSJ) = 



(_y+L+j+i Svi J6l(l + l)(2l + 1)(2A' + 1)(2A + 1)(2S> + 1)(2S + 1) 



X' 1 A 1 j S' 1 1 S \ 

X ( I L V )\ \ \ \ )\ X J X' ) 

x (P'L\ NL) (p'l'\ v(r) \nl) , (A10) 



(p'l' P'L X'S'J\ w(r)(o-i x <r 2 ) ■ I \nl NL XSJ) = 

(-) l+L+J+s 5 VI (5 s , 5 S i + Ss^Sso) v /l2Z(Z + l)(2/ + l)(2A' + l)(2A + l) 

(P'L\ NL)(p'l'\v(r)\nl) , (All) 



(A' 1 Al 


j S' 


1 




\ I L V ] 


I A 


J 


i) 



(p'l' P'L X'S'J\ v(r)r -V\nlNL XSJ) = 
(p'l' P'L X'S'J\ v(r)d s (r) \nl NL XSJ) = 

Si>i6 x . x 5s>s (P'L\ NL) (p'l'\ v(r)d s (r) \nl) , (A12) 

(p'l' P'L X'S'J\ v(r) <ti • a 2 r • V \nl NL XSJ) = 

(p'l' P'L X'S'J\ v(r) o-! • <t 2 d s (r) \nl NL XSJ) , (A13) 
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(p'l' P'L X'S'J\ v(r) a 1 ra 2 V \nl NL XSJ) = 

- (p'l' P'L X'S'J\ v(r) <r x ■ a 2 d s (r) \nl NL XSJ) 

+ - (p'l' P'L X'S'J\ v(r) S 12 (r) d T (l', I; r) \nl NL XSJ) 

+ - (p'l' P'L X'S'Jl i^- (o-i x a 2 ) ■ I \nl NL XSJ) , (A14) 
2 r 

(p'l' P'L X'S'J\ v(r) <r 2 f 01 • V \nl NL XSJ) = 

- (p'l' P'L X'S'J\ v(r) a x ■ a 2 d s (r) \nl NL XSJ) 

+ - (p'l' P'L X'S'J\ v(r) S X2 (r) d T (l', I; r) \nl NL XSJ) 
3 

- i (p'l' P'L X'S'J\ i-^- (<ti x a 2 ) ■ I \nl NL XSJ) . (A15) 
We have introduced the following effective radial differential operators: 

4M = J;, 

MO.) = | + ' (i + 1) -'J + 1) + 6 . (M6) 

When the short range correlations are implemented as indicated in Eqs. (J100I) . one should, 
accordingly, make the following replacements in the relative radial matrix elements in Eqs. ()A3|) - 

\nl) -» 5A^(^) l«0 , 

-» (pY| 5iV iv(r). (A17) 

Clearly, the correlation function gkN(r) is also subject to the action of the differential operators 
in these equations. Explicitly, the center-of-mass radial overlap and the relative radial matrix 
elements are given by 

(P'L\ NL) = J R 2 dRj L (P'R)1l NL (b/V2, Rj (A18) 

and 

(p'l'\v(r)d(r)\nl) = J r 2 dr j v (p'r) g NN (r)v(r) d(r) [ gAN (r) K n i (v^ b, r)] , (A19) 

where ji and jy are spherical Bessel functions, IZnl and 1Z n [ are harmonic oscillator radial 
wave-functions with the first arguments giving the respective length parameters, and d(r) is 
either unity or one of the differential operators (|A16|) . 

Finally, let us remark that the convention for the relative coordinate adopted here, Eq. (J14I) . 
has the opposite sign to those of Refs. |5J and jSJ, for instance. This introduces an extra phase 
factor of (—1)' and (—1)' , respectively, in our kets and bras involving the relative motion. As a 
result, the transformation brackets in Eq. (|A2|) differ by a factor of (—1)' from those originally 
defined in Ref. jSJ. To avoid this adjustment, one can simply shift to the opposite convention 
for the relative coordinate by making the transcriptions r — > — r, r — » —r and V — > —V in the 
expressions for the transition potential in coordinate space given in Subsections 12.31 and 12.41 In 
either case, the equations in Subsection 13. II and here in Appendix EI remain formally unaltered. 
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B Phase conventions for a, 6, . . . / 

With the phase conventions of Nabetani et al. (N.O.S.K.), in Ref. [HZ], from which Eq. (jl01|) is 
obtained, the transition amplitudes are defined as follows: 

a = ((pn ) 1 So|F|pA, 1 So)) NOSK , 
b = (<p«, 3 P |F|pA/So>) NOSK , 

d = ((pn, 3 D!|y|M, 3 Si)) 

V 1 1 VN.O.S.K. 

e = (<pn, 1 P 1 |^| P A, s S 1 >) N _ asjc , 

/ = ((p«, 3 p 1 |v r |M, 3 s 1 )) NOSK • (Bl) 

The relationship with our phase conventions is 

(<pn, ^I'jiripA, 25+1 ^)) NasK . = R 5 ' +5 ^' (np, 25 ' +1 6|^|Ap, 25 + 1 / J ) , (B2) 

where the first correction is due to the change in ordering in the Clebsch-Gordan couplings 
for the spins, and the second one, to the fact that we do not include the phase r in the final 
partial- wave radial function (cf. Eq. ()A19|) ). This explains the extra phase factors appearing in 
Eq. (QUI). 

Notice that there is no factor in Eq. jB^ as might be expected due to the change 

in ordering of the two particles in both the initial and final states, the reason being that this is 
compensated by the fact that we use the opposite convention for the relative coordinates. 
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